Valley dependent many-body effects in 2D semiconductors 
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We calculate the valley degeneracy (gv) dependence of the many-body renormalization of quasipar- 
ticle properties in multivalley 2D semiconductor structures due to the Coulomb interaction between 
the carriers. Quite unexpectedly, the dependence of many-body effects is nontrivial and non- 
generic, and depends qualitatively on the specific Fermi liquid property under consideration. While 
the interacting 2D compressibility manifests monotonically increasing many-body renormalization 
with increasing g^, the 2D spin susceptibility exhibits an interesting non-monotonic dependence 
with the susceptibility increasing (decreasing) with for smaller (larger) values of with the 
renormalization effect peaking around gv ~ 1 — 2. Our theoretical results provide a clear conceptual 
understanding of recent valley-dependent 2D susceptibility measurements in AlAs quantum wells. 
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Two dimensional (2D) electron (or hole) systems based 
on semiconductors, e.g. inversion layers, FETs, quan- 
tum wells, and heterostructures, have served as useful 
laboratory systems for studying electron-electron inter- 
action induced many-body effects for almost forty years 
The reason is the ease with which the 2D carrier 
density (n) can be tuned in these systems, thus effec- 
tively changing the 2D electron gas (2DEG) from being 
a weakly interacting system at large n to a strongly in- 
teracting system at low n. Such a variation in carrier 
density can not be achieved in 3D metals where only a 
modest tuning of the electron-electron interaction effect 
is possible by going from one metal to another with the 
concomitant complication of a varying background lat- 
tice structure [5]. The electron interaction strength in 
a many-body Coulomb quantum system is characterized 
Q by the dimensionless parameter which measures the 
average inter-particle separation in units of the effective 
Bohr radius. For the 2DEG, = {nn)-^/'^ / {nh^ /me^), 
where n,m, and n are respectively the 2D electron den- 
sity, the electron effective mass, and the background lat- 
tice dielectric constant. Often is considered to be the 
ratio of the average Coulomb potential energy, /{htq) 
where ro = (tth)"^/^ is the average inter-particle sepa- 
ration, to the average kinetic energy Ep = mrh? /m for 
a 2DEG (including the spin degeneracy). Although the 
appropriate definition of Vs, also called the Wigner-Seitz 
radius, is in terms of a dimensionless length, the defini- 
tion in terms of the dimensionless interaction energy is 
physically more appealing since > 1 (or < 1) region 
may be considered as the dilute strongly interacting (or 
the dense weakly interacting) regime. Varying the carrier 
density one can increase r^, thus accessing the strongly 
interacting dilute regime where the quasi-particle renor- 
malization correction due to many-body effects would be 
large. In practice, the very strongly correlated regime of 

« 10 ~ 30 [st] can be achieved in 2DEG whereas in 3D 
metals 1^], « 2 ^ 5. It is, therefore, not surprising that 



2D semiconductor systems have long served as the labo- 
ratory test-bed for studying r^ dependent interaction ef- 
fects. In general, the quasi-particle many-body renormal- 
ization is enhanced with increasing r^, and in principle, 
there could be quantum phase transitions (e.g Wigner 
crystallization, ferromagnetic instability, dispersion in- 
stability, etc.) in the strong-coupling Ts ^ 1 regime. 

Although electron interaction effects at T = are com- 
pletely characterized by the single dimensionless density 
parameter in a single- valley system (e.g. GaAs), a 
multi-valley semiconductor (e.g. Si, Ge, AlAs) is far more 
complex since the valley degeneracy (gv), i-e the num- 
ber of the equivalent valleys the electrons occupy in the 
ground state due to the semiconductor band structure, 
becomes an additional relevant parameter characterizing 
the electron-electron interaction strength. (We consider 
only equal valley population of all valleys and equal 
spin population of spin up-down levels at T = in this 
work, i.e. we consider only a paramagnetic ground state 
in both valley and spin quantum numbers.) It is obvious 
that in the presence of a valley degeneracy (i.e. an arbi- 
trary value of g.v > 1), both Ts and g^ will determine the 
many-body renormalization effects, but a careful inves- 
tigation of how gy itself affects the quantum many-body 
effects in 2DEG has not yet been carried out in the the- 
oretical literature [4] . This is precisely what we report in 
this Letter, concentrating on the interacting 2D spin sus- 
ceptibility and 2D compressibility and calculating their 
many-body renormalizations as functions of both and 

A relatively straightforward interpretation of the 
parameter as a dimensionless coupling energy, i.e, the 
Coulomb energy divided by the Fermi energy, gives = 
Coulomb energy = {e^ / Kr^) / {-Kh^n/ g^jin) = r^gv 
As stated before, rg = Ve for = 1, but for > 1, 
Tg > Ts, implying that the interpretation of the quantum 
Coulomb coupling as a dimensionless energy would imply 
that increasing g^ automatically involves a linear increase 



2 



of the dimensionless coupling parameter. In fact, a third 
possible definition of the dimensionless coupling strength 
is Tj. — y/g^rg , where now measures the dimensionless 
strength of the average exchange energy to the average 
kinetic energy at T = 0. Again, Tx increases with increas- 
ing gy at a fixed Ts implying that the quasiparticle renor- 
malization should be enhanced with g^. Although these 
alternative definitions of the Coulomb coupling constant 
(e.g. Tg, Tx) have often been emphasized in the literature 
[5| in order to claim that a multi-valley (i.e. g^ > 1) 
system is more strongly interacting than a single-valley 
system at the same total electron density (i.e. fixed r^), 
we assert that the many-body problem for the multi- 
valley situation is necessarily a two-parameter problem 
(i.e. Tg and g„) which can not be described in any situa- 
tion by a single effective parameter such as r,, = g^Vg or 

= ^Jg^Tg or any other combination of Vg and g^. The 
non-trivial dependence of the interaction-induced many- 
body effects on the two independent parameters Vg and 
gv-i in particular, how even the qualitative nature of the 
quasiparticle renormalization as a function of Vg and g^ 
manifests itself completely differently in different proper- 
ties of the multivalley 2DEG (e.g compressibility versus 
susceptibility), is the central theme of this work. Specif- 
ically, we show that many-body effects could either be 
enhanced or suppressed with increasing g^ . 

Although the main motivation of our work is theoreti- 
cal, a part of our motivation comes from the extensive re- 
cent experimental work on the multivalley AlAs 2D quan- 
tum wells carried out at Princeton university 0]. These 
experiments demonstrate that the measured 2D spin sus- 
ceptibility of AlAs quantum wells depends on g^ (as well 
as Ts), and in general the susceptibility is smaller for 
larger values of g^ with the difference between g^ = 1 and 
2 decreasing with decreasing Vg (i.e. increasing density). 
It was emphasized in these experimental papers ^ that 
such a higher value of susceptibility for gy = 1 compared 
with ^i, = 2 is contrary to the popular wisdom, based on 
considerations which claim r^i— gv^s) or rx{— ^/g^fg) to 
be the appropriate interaction parameter, which would 
imply an increasing many-body renormalization with in- 
creasing gy. In the current work, we resolve this puzzle 
by showing that the many-body enhancement of the 2D 
spin susceptibility decreases with increasing gy (between 
I and 2) except at very small values of rg. We also make 
the predication that an increasing g„ will always increase 
the many-body renormalization effect of the compress- 
ibility (for all values of Vg and gy) in a sharp contrast 
with the susceptibility. 

We employ the one-loop self-energy calculation as our 
basic underlying theory to calculate the T — suscepti- 
bility and compressibility as a function of rg and gy . The 
self energy is being calculated in the leading order ex- 
pansion in the dynamically screened Coulomb interaction 
because of the well-known long distance (i.e. q ^ 0) di- 
vergence of the bare Coulomb interaction. This is equiva- 



lent to calculating the thermodynamic grand potential in 
the infinite ring or bubble diagram expansion and then 
obtaining the spin susceptibility and the compressibil- 
ity by taking the appropriate derivatives with respect 
to the magnetic field and the volume respectively. The 
important thing to remember is that each bubble dia- 
gram carries a factor of gy due to the valley degeneracy 
(in addition to the factor of 2 due to the spin degener- 
acy), and the Fermi wave vector kp and the Fermi energy 
Ep (i.e. the non-interacting chemical potential) are both 
suppressed by the valley degeneracy: 



kp — (2Trn/gyy^^;Ep = Trh'^n/{mgy) 



(1) 



where n is the 2D carrier density and m is the bare ef- 
fective mass. The basic 2D non-interacting polarizibility 
function, i.e. the bubble diagram, is given by, 

U{q,cu) = -2gy f (2) 

J (27r)^ ftw - [ek+q - Sk] + i5 

where Sk = h^k^ /2m, and fk is the Fermi distribution 
function corresponding to the energy E{k). We note that 
the 9 ^ 0, w ^ limit of the bubble diagram gives the 
2D non-interacting density of states, D{E) = gym/nh^ , 
which is enhanced by > 1. Eqs. ([T]) and ([2|) show how 
gy enters the theory: (I) By decreasing Ep (Eq. ((T|)), 
gy tends to enhance interaction effects; (2) by increas- 
ing screening through Eq.(l2]), i.e. by enhancing the den- 
sity of states, gy tends to suppress the interaction ef- 
fects. Depending on which of these effects is more im- 
portant in determining a particular property, the many- 
body renormalization may increase, decrease, or show a 
non-monotonic behavior with increasing gy. 

We skip the technical details of the single-loop self- 
energy calculation, and provide below the final formula 
we use for calculating the interacting susceptibility (x) 
and compressibility {K) in terms of rg and gy, 

dq ^ 
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du- 



e{q,iuq) 



dk + b'^ 



k=l 



Ko 
K 



1 d{kl/2 + Re[i:{kp,Ep)]) 



(3) 



(4) 



kp dkp 

In Eqs. ^ and (H]), xo and Kq refer to the corresponding 
non-interacting quantities, a = + /c^ — 9^/4, b = uq, 
and e{q, iuq) = 1 — Vqll{q, iuq) is the dynamical RPA di- 
electric function, and T,{kp,Ep) is the one-loop self en- 
ergy calculated at the Fermi surface, and a = ^ gv/2 and 
kp = ll{arg). We note that the final equations are sufh- 
ciently complex that it is not possible to read off the {rg , 
gy) dependence for arbitrary values of rg and gy. Clearly, 
the results depend on both rg and gy independently. We 
have, therefore, numerically calculated the (rg, gy) de- 
pendence of the susceptibility and the compressibility for 
arbitrary rg and gy values. 
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FIG. 1; (a) Susceptibility vs Ts for different valley degeneracy. 
The solid line is for the numerical data, dash-dotted line for 
the asymptotic formula shown as Eq. ([5]). Inset gives the in- 
verse spin susceptibility vs Ts < 1 showing non-monotonicity. 
(b) Susceptibility vs > 1 for different valley degeneracy. 



First, we present our susceptibility results in Fig. 1 
and 2. In Fig. 1, we show xlxo ^-s ^ function of (for 
— 1,2,4,6) in the small Ts regime (r^ = — 0.2) 
in Fig. 1(a) and in the large Vg regime {ts — \ — 7) in 
Fig. 1(b). It is clear that the situation is qualitatively dif- 
ferent between the small (r^ < 1) and the large (r^ > 1) 
rg regime. For small Vg, the many-body effects are en- 
hanced by increasing in accordance with one's naive 
expectation based on the suppression of Ep by g^. For 
large values, on the other hand, the many-body renor- 
malization of x decreases with increasing g^, in complete 
agreement with the recent experimental measurements 
[6| which were carried out in the > 1 regime. One 
of our predications is, therefore, that a measurement of 
the 2D susceptibility at high densities (r^ < 1) would 
manifest a non-monotonic crossover in the gv dependence 
of the many-body effects as can be seen in the inset of 
Fig. 1(a). In Fig. 1(a), we show by the dotted Hue our an- 
alytic formula valid in the <C 1 regime which can be 



FIG. 2: (a) Susceptibility vs for different electron density 
Ts. For small values, the non-monotonicity can be seen in the 



inset, (b) Inverse susceptibility vs I/Qv for different rs 
shows the constant x for <7jj 3> 1. 
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derived to be 



9v 



2.22 



rjgj 

2V27r' 



(5) 



We can see from Eq. ([5]) that x can not be written, even 
in the <C 1 limit, as a single-parameter function. 

In Fig. 2, we show the g^ dependence of the interact- 
ing susceptibility, both for small and large g^. The non- 
monotonic behavior of the susceptibility as a function 
of gv for smaller values is obvious in Fig. 2, and the 
decrease of x with increasing g^, eventually reaching a 
constant for unphysically large g^ (all g„ > 6 is unphysi- 
cal since no known multivalley 2D semiconductor system 
has more than six valleys). For < \, the maximum 
in x/xo occurs between gv — ^ and 2, and should there- 
fore be experimentally accessible as a matter of principle. 
For Ts > 1, however, the maximum in Fig. 2 moves to 
the gv < 1 regime which is unphysical. The very large 
ffi; ^ 6 regime (Fig. 2b) is of theoretical interest only. 
We obtain the following asymptotic analytic formula for 
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FIG. 3: (a) Compressibility vs Vs for diflerent valley degener- 
acy, (b) Compressibility vs (?„. Dotted lines are HF results. 
Inset shows the compressibility vs for large gv 
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(6) 



up to logarithmic corrections in 
formula (Eq. ([6])) agrees well with our numerical results 
for ^ 100. The small (and the small Ts) regime 
shown in Fig. 2 agrees with the analytic Eq. ([5|) derived 
for the small Ts regime, where the non-monotonicity with 
respect to can be seen to be arising from the negative 
sign in the last term. We see that Eq. ([5]) implies a max- 
imum in X at 5™ ~ 1.07 (for = 1) and 1.78 (r^ = 0.5), 
which is consistent with the numerical result of Fig. 2. 

In Fig. 3, we show our calculated Ts and g^ depen- 
dence of the 2D interacting (inverse) compressibility. It 
is clear that, in contrast to the interacting susceptibility, 
the interacting compressibility manifests monotonically 
stronger many-body effects with increasing valley degen- 
eracy g^ in all regimes of Ts and g^. There is, in fact, no 
non-monotonicity at all in the interacting compressibil- 
ity which decrease continuously with increasing or , 



eventually becoming negative for Ts ^ 2.5 in a 2DEG 
We can analytically calculate the small dependence of 
the interacting compressibility, obtaining: 



Kirs < 1) 



1 - 0.45r,gy2 _^ o,022{r%) logir^gl^^), 

(7) 

where we note that the first two terms of the expan- 
sion (i.e. the Hartree-Fock result with just the exchange 
self-energy correction) are identical to the corresponding 
first two terms (i.e. the Hartree-Fock result) for xo/x i^i 
Eq. ([5|). We mention that the leading-order correlation 
correction, the third term in Eq. has an r^gv depen- 
dence in the compressibility (in contrast to the r^g^ de- 
pendence in Eq. ([5]) for the susceptibility) . Consequently 
the correlation correction to the compressibility is very 
weak for small values. In Fig. 3(b), we compare 
the Hartree-Fock result with our numerical results, and 
for small r^, the agreement is excellent. For large g^, 
Kq / K decreases monotonically linearly in g^ (neglecting 
logarithmic corrections) . 

In discussing the significance of our results, we empha- 
size that we have provided a complete qualitative reso- 
lution of the puzzle raised by the experimental observa- 
tions of ref. [6|], where the 2D interacting susceptibility 
was found to be smaller for the larger valley degeneracy 
system. In addition, we make a clear predication that 
if the experiments are carried out in higher-density (i.e. 
lower Ts) samples than used so far, the susceptibility will 
be larger in the higher valley degeneracy system. We also 
predict that the many-body effects in the 2D compress- 
ibility, in sharp contrast to the 2D susceptibility, would 
be enhanced monotonically with increasing g^, for all 
values. We mention that in our theory the valley (g^) and 
the spin (g^) degeneracy are equivalent, and therefore a 
measurement of valley susceptibility in the presence of 
a variable spin degeneracy (e.g. gs = 2 or 1) would mani- 
fest exactly the same qualitative trend as seen in Figs. 1 
and 2 with the roles of valley and spin being interchanged. 
Similarly we predict a suppression of the compressibility 
with decreasing spin degeneracy from ga — "2, to Ihy ap- 
plying an external parallel magnetic field at fixed value 
of g^, ■ We believe that an experimental confirmation of 
our qualitative predictions would be compelling evidence 
of the Fermi liquid nature of 2D interacting systems even 
at large Ts and g^ values. 
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